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Abstract. The interplay between dust and molecules is of fun- 
damental importance in early galaxy evolution. First we present 
the prediction for the dust emission from forming galaxies. Then 
we discuss the observational strategy for molecules in early galax- 
ies by infrared absorption lines of a bright continuum source be- 
ly-j hind the clouds. By combining these two approaches, we will be 

£f~) | able to have a coherent picture of the very early stage of galaxy 

00 ■ evolution. 

On 

O . 

*0 ■ 1 Introduction 

O ; 

C"| . At the very first phase of the Universe, there have been almost no struc- 
i' tures, as seen in the tiny fluctuation of the cosmic microwave background. 
q , In contrast, today we observe a variety of cosmic structures, from planes 
!— i ' to the large-scale structure in the Universe. In addition, the Big Bang nu- 
c/2 , cleosynthesis produced only light elements, while the present Universe is 
rich in heavy elements. Thus, galaxy formation may be the most important 
event in the cosmic history. 

Active star formation is followed by heavy clement production from the 
birth and death of stars. The produced heavy elements generally exist in 
' the form of small solid particles, i.e., dust. Dust grains play very important 
roles in galaxy evolution. First, they scatter and absorb short wavelength 
photons and finally re-emit the energy in the infrared (IR). This makes 
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the galaxy evolutions seen in the optical/ultraviolet and IR very different 
with each other (see Takeuchi et al. 2005c) . Another important role of dust 
grains is that they work as a catalyst of molecular formation (e.g., Hirashita 
et al. 2002, Hirashita & Ferrara 2002). This process is particularly impor- 
tant in the very early phase of galaxy evolution, because hydrogen molecules 
are very important coolant of gas in the very metal-poor phase of galax- 
ies. Without dust, star formation does not proceed effectively (Hirashita & 
Ferrara 2002). 

From this point of view, we focus on the two kinds of systems. We first 
discuss young galaxies with active dust production. We will see that their 
appearance is different from dusty IR luminous galaxies at lower redshift 
(e.g., Takeuchi et al., 2001a, b). We discuss the observability of their con- 
tinuum radiation from dust. Then we move the topic to dense gas systems 
with little metal/dust. On the eve of the active dust production, they might 
contain a large amount of hydrogen molecules. We consider the direct mea- 
surement of the hydrogen molecules in such protogalactic clouds through 
the absorption lines in the IR. 

2 Dust Emission from Forming Galaxies 

2.1 SED Model for Forming Galaxies 

We present a brief outline of our model framework. All the details are shown 
in Takeuchi et al. (2003, 2005b), and Takeuchi & Ishii (2004). 

Nozawa et al. (2003) investigated the formation of dust grains in the 
ejecta of Population III SNe, whose progenitors arc initially metal-free. 
They considered unmixed and uniformly mixed cases in the He core. In 
the unmixed case, the original onion-like structure of elements is preserved, 
and in the mixed case, all the elements are uniformly mixed in the helium 
core. 

For the chemical evolution, we use a closed-box model with the Salpeter 
initial mass function 0(m) oc to~ 2,35 with mass range of (mi, m u ) = (0.1 M , 100 
We neglect the contribution of SNe la and winds from low-mass evolved stars 
to the formation of dust, because we consider the timescale younger than 
10 9 yr. The interstellar medium is treated as one zone, and the growth of 
dust grains by accretion is neglected. Within the short timescale considered 
here, it can be assumed safely. We also neglect the destruction of dust grains 
within the young age considered. Finally, we assumed a constant SFR for 
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simplicity. 

Very small grains cannot establish thermal equilibrium with the ambient 
radiation field, which is called stochastic heating. For the specific heat, 
we adopt a multidimensional Debye model (e.g., Draine & Li 2001) for 
carbon and silicate grains. For other species, we adopt the classical three- 
dimensional Debye model with a single Debye temperature. The emission 
from dust is calculated basically according to Draine & Anderson (1985). 
Total dust emission is obtained as a superposition of the emission from each 
grain species. We constructed Q(a, A) of each grain species from available 
experimental data via Mie theory. When the dust opacity becomes large, we 
should consider the self- absorption of IR emission by dust. This is treated 
by a thin shell approximation. 

2.2 Results and Discussions 

2.2.1 Evolution of IR SED 

We first show the evolution of the IR SED of forming galaxies based on 
our baseline model in Figure ^ F°r these calculation wc adopted the star 
formation rate SFR = lM Q yr _1 . We adopt tsf = 30 pc and 100 pc. In 
a very young phase (age = 10 6 ' 75 -10 7 25 yr), unmixed-case SED has an 
enhanced N-MIR continuum. After 10 725 yr, the N-MIR continuum is 
extinguished by the self-absorption in the case of tsf = 30 pc. In contrast, 
the self-absorption is not significant for rgp = 100 pc. In both cases, the 
SEDs have their peaks at a wavelength A ~ 20-30 /im. 

2.2.2 A nearby forming galaxy SBS 0335-052 

Although it is still difficult to observe galaxies in their very first phase of 
the SF, we can compare the model with a local counterpart of such galaxies. 
SBS 0335-052 is a local galaxy (~ 54 Mpc) with SFR = 1.7 Mgyr" 1 and 
extremely low metallicity Z = 1/41 Zq. This galaxy is known to have an 
unusual IR SED and strong flux at N-MIR. It has a very young starburst 
(age < 5 Myr) without significant underlying old stellar population. Houck 
et al. (2004) presented new data of the MIR SED by Spitzer. We have 
calculated the SED for r SF = 10, 20, and 30 pc. The SFR is fixed to be 
1.7 Mg yr -1 , and the age is 10 6 5 yr. We present the model and the observed 
SED of SBS 0335-052 in Figure The very strong N-MIR continuum of 
SBS 0335—052 is well reproduced by the SED of unmixed case, but the 
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Fig. 1. The evolution of the IR SED of a very young galaxy (left six panels: 
rsF = 30 pc; right six panels: 100 pc). 



mixed case seriously underpredicts its MIR continuum. This suggests that 
we may determine the dust production (unmixed or mixed) of SNc through 
the observation of the N-MIR SEDs of forming galaxies. 



2.2.3 Toward higher redshifts 

Consider a dark halo of mass ~ 10 9 M©, then it is expected to contain a gas 
with mass ~ 10 8 Mq. We calculate the SEDs for a dwarf forming galaxy. If 
gas collapses on the free-fall timescale with an efficiency of esF (we assume 
£sf = 0.1), we obtain the following evaluation of the SFR (Hirashita & Hunt 
2004): 
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Wavelength A [yLim] Wavelength A [/^m] 

Fig. 2. The model for the SED of a nearby star- forming dwarf galaxy 
SBS 0335-052. 

If we consider M gas ~ 10 s M®, we have SFR ~ 3(r SF /100 pc)~ 3 / 2 Mgyr" 1 . 
Thus, we consider a dwarf galaxy with SFR = 10 M Q yr _1 as an example, 
and we adopt r$F — 30 pc and 100 pc. The age is set to be 10 7 yr. We 
show the expected SEDs for such galaxies at z = 5, 10, and 20 in Figure 
It seems almost impossible to detect such objects by Herschel or ALMA. 
However, gravitational lensing works very well as a natural huge telescope. 
This is depicted by the thin lines in Figure|3|(with magnification factor 40). 
If we suppose a cluster of galaxies at z\ ~ 0.1-0.2 whose dynamical mass 
■Mdyn is 5 x 10 14 Mq, and use number counts of high-z forming galaxies 
presented by Hirashita & Ferrara (2002), we have an expected number of 
galaxies suffering a strong lensing to be ~ 1-5 with detection limit 1 fjjj. 

3 Direct Measurement of H2 by IR Absorption Lines 

3.1 Basic Idea 

Molecular hydrogen is the predominant constituent of the dense gas in the 
Universe. While in the local Universe, molecules containing heavy elements 
(e.g. CO, H 2 0) are good tracers of the amount of molecular hydrogen, they 
are expected to be significantly depleted in the early Universe. We need 
a technique for measuring molecular hydrogen directly. Petitjean et al. 
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Fig. 3. The expected flux densities of a dwarf forming galaxy at z = 5, 10, and 20. 
Gray and black tick marks are the confusion limits of Herschel and the detection 
limits of an ALMA 8-hour survey. The Herschel confusion limits are taken from 
Lagache et al. (2003). 



(2000) have reported a direct measurement of H2 molecules in the UV. The 
transition probabilities of ionization and dissociation lines arc so large that 
they are useful for detecting thin layers and small amounts of the molecular 
gas, but useless for detecting dense gas clouds. 

In contrast, molecular hydrogen has well-known vibrational and rota- 
tional transitions at IR wavelengths. Their transition probabilities are very 
small because hydrogen molecule, a diatomic molecule of two identical nu- 
clei, has no allowed dipole transitions. Hence, they are useful tools to an- 
alyze dense (> 10 cm~ 3 ) and hot (> 300 K) gas. Ciardi & Ferrara (2001) 
calculated expected emission line intensities of molecular hydrogen, but un- 
fortunately, their direct measurement was found to be very difficult due to 
their weakness. However, if there is a strong IR continuum source behind 
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or in the molecular gas cloud, absorption measurements of these lines may 
be possible. This idea is presented by Shibai et al. (2001). 

We should note that the absorption by dust in the cloud is usually 
larger than the absorption by H2 molecules. However, since our target 
is a primordial gas cloud whose metallicity is significantly lower than the 
Galactic value, the lines can be detected in absorption against bright IR 
sources. Such observation will be feasible with the advent of proposed space 
missions for large IR telescope facilities. Here we focus on SPICA (Ueno et 
al. 2000). 

3.2 Calculation 

Assuming a uniform, cool gas cloud (kT CK <C hv), the optical thickness of 
the line absorption, rii nc is 



where the subscripts u and I indicate the upper and lower levels of a tran- 
sition, g u and gg are the degeneracy of each state respectively, A u g is the 
Einstein's A coefficient, Ni is the column density of the molecules in the 
lower state, and AV is the line width in units of velocity. Here we as- 
sume that almost all the molecules occupy the lowest energy state. The 
absorption line flux in the extinction free case, Ifi^l o> is obtained by 



where Av is the line width in units of frequency and S is the continuum 
flux density of the IR source behind the cloud. Assumed parameters in 
this calculation are listed in Table ^ Next, we consider about the dust 
extinction which is denoted by 




(3.1) 




(3.2) 




(3.3) 



The extinction spectrum of Mathis (1990) is adopted for (Ax/Ay), and the 
extinction efficiency is assumed to be proportional to the relative heavy 
element abundance, and (Av/N-r)q is the conversion factor from Ay to 
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Table 1. Parameters assumed for the present calculation 

Ortho : Para 3 : 1 

HD/H 2 1CT 5 

Intrinsic Flux of Source 10 mJy 

Line width 100 km s -1 

Detection Limit (SPICA 5cr) 5 x 10~ 21 Wm~ 2 

Detectable Optical Thickness > 0.01 

Dust Extinction Model Mathis (1990) 



Nh for the local abundance. We then obtain the absorption line flux with 
extinction, Jgjjjj, 

I&l = Cl' ^^ = SA » (1 - e " Tlinc ) e_TdUBt ■ ( 3 - 4 ) 



3.3 Results 




Fig. 4. Left panel: the absorption line flux expected for a cloud with the Milky 
Way heavy element abundance (Z = 1) in front of a 10 mJy source. Right panel: 
same as (a) but the heavy element abundance is Z = 0.01. 

The calculation was made for the five lines listed in Table [21 Left panel 
of Figure 01 shows the absorption line flux expected for a cloud with the 
local heavy element abundance in front of a f mJy source. The absorption 
fluxes are far smaller than the 5er expected sensitivity of the SPICA mission 
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Table 2. Parameters of the Line Transitions 



Transition Wavelength A coefficient A\j Ay 



H 2 v = 


1 - 


5(1) 


2.12 


3.47 x 10- 


-7 


0.11 


H 2 v = 


1 - 


5(0) 


2.22 


2.53 x 10" 


-7 


0.11 


B 2 v = 


- 


5(1) 


17 


4.77 x 10- 


-10 


0.020 


H 2 v = 


- 


5(0) 


28 


2.95 x 10" 


-11 


0.011 


HD v = 


0- 


R(0) 


112 


2.54 x 10" 


-8 


0.0011 



(Ueno et al. 2000). On the other hand, right panel of Figure 0] shows the 
result for the case in which the heavy element abundance is 1 % of that 
of the local one. All five lines populate parameter space above or near the 
limit. 

We have assumed an intrinsic continuum flux density of 10 mJy at the 
line wavelength. We examined some existing QSOs as template sources (for 
the details, see Shibai et al. 2001). If we have objects at z — 5 at least as 
luminous as our template sources, it is possible to observe absorption in H 2 
lines of 17, 28, and 112 /im respectively. The redshifted 112 /in line will be 
in the submm range and may not be observable by IR missions, but ALMA 
should be able to observe it. Unfortunately, it is very difficult to observe 
the absorption of 2.2 /im lines since the expected flux density based on dust 
emission spectrum will be too faint. 



3.4 What can we learn from the absorption lines of these clouds? 

First we simply estimate the column density of a protogalactic hydrogen 
cloud. Here we assume that the size of the cloud R is a few kpc and that 
the gas mass M is ~ 10 11 M Q . Such a large reservoir of molecular gas has 
been discovered at high redshift (Papadopoulos et al. 2001). Such clouds 
have column density of 

where / is the mass fraction of the molecular clouds to the total gas mass, 
and M is the total gas mass of the protogalaxy. The radius will be finally a 
few kpc. In a realistic situation, a primordial cloud may evolve dynamically 
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on a free-fall timescale, which is much shorter than the timescales of cos- 
mological structure evolution. Therefore, observed properties are specific to 
the redshift at which the cloud absorption is measured. Hence, since we will 
obtain redshift z, and velocity dispersion AV, we can trace the dynamical 
evolution of a primordial cloud at high-z. The SPICA mission will provide 
us with the information on the massive objects (> 10 11 M Q ) at z < 5. We 
still have to wait for more sensitive facilities for the observational approach 
to the Population III objects, whose physical properties are theoretically 
predicted by recent extensive investigations, because they are located at 
higher redshift and their typical mass is small (e.g., Nishi & Susa 1999). 
When wc estimate the detectability of such objects by other facilities, the 
algorithm constructed in this paper can be applied straightforwardly. 

4 Conclusion 

We first present our dust emission model from forming galaxies, which is 
based on a new theory of SN II dust production. The model roughly re- 
produced the observed SED of a local low-metallicity dwarf SBS 0335—052 
which has a peculiar strong and MIR-bright dust SED. We also calculated 
the SED of a very high- z forming small galaxy. Although it may be intrin- 
sically too faint to be detect ed even by ALMA 8-hour survey, the gravita- 
tional lensing can make it feasible. 

Then we proposed a method to measure the amount of H2 in primordial 
low-metallicity cloud in an IR spectra of QSOs. If metallicity of the cloud 
is low (Z ~ 0.01Z Q ), dust extinction is so weak that 17 and 28 /xm lines are 
detectable by SPICA for objects at z < 5. For very high-z Population III 
objects, ALMA will be useful. By this observation, we can trace back the 
dynamical evolution of early collapsing objects at very high-z. 

By combining these two approaches, we will be able to have a coherent 
picture of the very early stage of galaxy evolution, from the pristine gas 
phase to active dust production phase. 

This work is made through collaborations with many people, especially 
Andrea Ferrara, Hiroyuki Hirashita, Leslie K. Hunt, Takako T. Ishii, Takashi 
Kozasa, Takaya Nozawa, and Hiroshi Shibai (alphabetic order). I have 
been supported by the Japan Society of the Promotion of Science as a 
Postdoctoral Fellow for Research Abroad (Apr. 2004-Dec. 2005). 
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